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A  POINT  EXPLOSION  IN  AN  INHOMOGENEOUS  ATMOSPHERE 

^/Following  is  a  translation  of  an  article  by  A.  S.  Komparieneyets 
in  the  Russian-language  periodical  Doklady , Akademii  nauk  SSSR 
(Reports  of  the  Academy  of  Sciences  USSR),  Moscow,  VolT  130, 

No.  5,  pages  1001-1003^ 

As  is  well  known,  the  problem  of  a  point  explosion  in  a  medium 
of  constant  isentropic  ,.ndez  y  ^as  an  analytical  solution^ )  for 
types  of  phenomena  in  which  it  is  possible  to  neglect  the  energy  in 
the  medium  up  to  the  time  of  the  explosion  with  respect  to  the  energy 
carried  in  the  detonation  wave*  A  relatively  simple  solution  can  be 
obtained,  because  the  problem  can  be  put  into  a  dimensionless  forms 
it  contains  no  liharact eristic  parameters  with  the  dimensions  of  length*, 
velocity,  or  time.  If  there  is  a  density  gradient  in  the  medium,  even 
if  only  in  one  direction, it  is  dimensionless  formulation  is  no  longer 
possible  and  an  exact  solution  is  not  obtained.  For  small  gradients 
the  problem  can  be  considered  using  a  pertubation  methcd(2,3). 

If  the  explosion  occurs  in  an  extremely  rarefied  atmosphere, 
a  strong  shock  wave  with  a  limiting  compression  is  propagated  for 
such  a  distance  that  the  density  is  changed  by  many  times  with  respect 
to  the  density  at  the  site  of  the  detonation.  In  this  case  it  is 
insufficient  to  lineraze  the  equations*  The  numerical  solution  of 
the  problem  with  three  variables  is  very  laborious,  even  using:  an 
electronic  computer, 

A  s  emiqual  it  a  t  iv  e  method  may  bo  proposed  that  is  based  on  ;aii 
essential  property  of  the  exact  spheric all v  symmetric  solution.  In 
this  solution  the  energy  is  found  to  be  propagated  almost  uniformly 
through  the  whole  volume  of  the  detonation  wave,  and  it  is  only  close 
to  the  front  itself  that  it  exceeds  by  a  factor  of  2  or  3  the  mean 
value  throughout  the  volume.  The  whole  mass  of  the  material  is  also 
distributed  concentrically  throughout  ■'this  region. 

It  is  natural  to  assume  that  a  detonation  wave  in  an  inhomo-  , 
geneous  atmosphere  also  possesses  this  property.  Actually,  if  the 
pressure  within  the  wave  is  constant  in  space  (the  pressure  is  pro¬ 
portional  to  the  energy  density),  and  the  mass  density  equals  zero, 
the  hydrodynamic  equations  in  the  main -volume  are  trivially  satisfied. 
Thus,  in  order  to  describe  the  propagation  of  the  wave,  it 1  is  necessary 
vo  use  the  conditions  in  the  shock  front  itself. 

If  the  equation  of  the  wave  front  is  f(r,z,t)  •  0,  the  norfcal* 
vol 00 ity  component. -of  the*  front  Dn  is  defined  by  the  well  known 
equation,  -  :  : 
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Hero,  els  is  usual  in  tho  consideration  of  strong  waves,  the  initial 
pressure  is  neglcctod  with  respect  to  tho  pressure  of  tho  wave  front 
p.  In  this  approximation  the  density  behind  the  front  p*  is  con¬ 
nected  with  the  density  before  the  front  p  by  the  constant  relation¬ 
ship 


The  pressure  is  expressed  in  terms  of  the  energy  density  C 

p  s  (Y-l)€=  (V  —  (3) 

V 


whore  E  is  the  energy  density  of  the  detonation;  V  is  the  volume 
occupied  by  the  detonation  wave;  and  a  coefficient 

indicating  how  many  times  larger  the  energy  density  near  the 
front  is  than  the  mean  density  throughout  the  volume  (•*•).  The 
assumption  that  A  is  constant  over  the  surface  is  the  basis  of  the 
method  hfero  proposed. 

We  shall  consider  that  the  equation  for  the  wave  front  in 
cylindrical  coordinates  is  solved  for  tho  radius :  r  «  r(z,t).  Then 
the  total  volume  of  the  wave  is 

V(t)  =  n  [  r2  (z,  t)  dz,  (4) 

V 

whore  r(zi,  t)  »  Z£,  t)  =  0,  Substituting  (2),  (S),  and  (4)  in  (l), 
and  expressing  tho  density  by  tho  barometric  formula,  we  obtain  a 
partial  differential  equation  for  tho  function  r: 


(21  V 
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*  0.  (5) 


Hero  z0  is  the  equivalent  thickness  of  the  atmosphere,  and  y  is  an 
auxiliary  variable  defined  by  the  oquation 

t 
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P0  is  the  initial  air  density  at  the  site  of  the  detonation  (z  =  0)# 
The  equation  (5)  is  soluble  by  the  method  of  separation  of 
variables 

s  _ _ 

r  ■  fy  +  (  W?  *•  “  *Ao  —  l!  (7) 
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z 

’■?  r  L  h  e —  z/z° 

—  -  y  -f-  \ ,  d«  r=- «  r  (.- ),  (s) 

o  Wf  e  z/z0  — -  1 

For^ smo.ll  t  or  y  the  wave  may  be  spherical.  For  this  to  happen 
it  is  sufficient  to  put  the  function  F(;)  equal  to  zero.  Then  by 
eliminating  ”  from  (8)  and  substituting  in  (7),  we  got 


j  i 

r  =  2z0  arc  cos  !  _1  Gz/2z0  (  i  —  x2  ,  G_  z/z0)  \.  (9) 

i _ 2  _J 

hero  x  =  y/2zQ. 

Thus  we  obtain  the  position  of  tho  upper  and  lower  points  of 
the  wavo  z-^  and  Zg 

e  *l-»2/2z0  s  rp  (io) 

as  well  as  the  position  and  magnitude  of  its  maximum  radius 


0—  Vzo  ”  1  r—  rm  —  2z0  arc  sin  x.  (11) 

Thus  tho  maximum  possible  radius  of  the  wave  is  oqual  in 
this  case  to  x  s  ls  so  that  tho  upper  limit  of  thowwavo  goes  out  to 
..•nfinity*  But  this  occurs  after  a  finite  time  ^  which  is  dofinod 
from  (6)  as 

- -  i 

8l,Z0Po 
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-2 {a(l*x)  '  v 

■  J\  (x)  \  du  arc  cos2  j  _L  oV2  (i  «  x2  >(.0-u)  j.  (13) 

-2  In(l  +  x) 


The  timo  for  the  upward  wave  to  go  to  infinity  remains  finite, 
thanks  to  the  fact  that  the  volocity  of  tho  wavo,  according  to  (l), 
tends  to  infinity  as  z - *s  00.  Fig.  1  shows  the  curves  for  the  trans¬ 
formation  from  x  to  t  and  tl/2.  Drawn  to  soalo  in  Fig.  2 

are  the  calculated  cross  sections  of  a  wavo  in  a  vertical  plane 
passing  through  tho  site  of  tho  detonation,  for  various  instants  of 
t  ime . 

It  is  obvious  that  tho  solution  loses  its  significance  boforo 
z\  becomes  infinite,  Nonetheless  it  is  possible  to  rnako  the  follow¬ 
ing  conclusions  however  large  the  total  cnorgy  of  tho  detonation, 
a  strong  shock  wave  can  be  propagated  downward,  according  to  tho 
law  obtained  here,  by  not  moro  than  1.38  Zq  or  about  11  km.  With 
further  downward  propatation,  tho  shock  front  will  bo  attenuated 
moro  quickly  than  tho  following  wave  on  account  of  tho  rarefaction 
waves  proceeding  upward  from  it  into  empty  space.  Tho  propagation 
of  a  wave  through  unperturbed  air  will  recall  the  short  shock  in  a 
substance  adjoining  a  vacuum,  which  was  considered  by  Ya.  B*  ZolT 
dovich  (4). 
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